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We theoretically analyzed a magnetic wall confined in a nano-contact spin valve paying special 
attention to the penetration of the magnetic wall into thin ferromagnetic electrodes. We showed 
that, compared with the Bloch wall, the penetration of the Neel wall is suppressed by increases 
of the demagnetization energy. We found the optimal conditions of the radius and height of the 
nano-contact to maximize the power of the current-induced oscillation of the magnetic wall. We 
also found that the thermal stability of the Bloch wall increases when the nano-contact's radius 
increases or height decreases. 



Magnetic walls have attracted much attention as a ba- 
sic element of nano-spintronics devices such as a race- 
track memory^, spin-wave logic gates^, a read-head for 
ultra-high-density magnetic recording 3 - - — , and a spin- 
torque oscillator (STO)£~— . The nano-contact spin valve 
with a geometrically confined magnetic wall shown in 
Fig. Q] (a) is a nanostructure suitable for the latter two 
applications, and much effort has been devoted to study- 
ing its electric and magnetic properties'^ - — . It is 
widely accepted that a magnetic wall should be created 
in the nano-contact when the magnetization vectors of 
the top and bottom electrodes arc aligned to be anti- 
parallel. The magnetic wall in the nano-contact can be 
categorized into two types: the Bloch wall and the Neel 
wall. The magnetic moments of the Bloch (Neel) wall ro- 
tate along the axis parallel (perpendicular) to the nano- 
contacfeil~— . The read head utilizes the resistance change 
due to the creation of a magnetic wall, and the STO uti- 
lizes the oscillation between the Bloch wall and the Neel 
wall induced by the applied direct current. 

In early experiments 3 - ' 12 ' 16 the break junction of a fer- 
romagnetic wire was employed to realize the nano-contact 
spin valve, and the corresponding theoretical analysis 
predicted that the magnetic walls were almost completely 
confined in the contact region. The point is that the 
energy required to rotate the magnetic moments in the 
contact is much less than that in the electrode of the 
ferromagnetic wire. However, if the thickness of the elec- 
trode is as small as the height of the contact, the energy 
required to rotate the magnetic moments in the electrode 
becomes so small that the magnetic wall can penetrate 
into the electrode^— as shown in Fig Q] (b) . Recently 
the nano-contact spin valve with thin ferromagnetic elec- 
trodes^ 3 - has been considered a powerful candidate for 2-5 
Tb/in 2 read sensors 2 ^. Therefore, it is important to an- 
alyze the penetration of the magnetic wall into the thin 
ferromagnetic electrodes of the nano-contact spin valve. 

In this paper, we derived analytical expression of a 
magnetic wall penetrating into the thin ferromagnetic 
electrodes of a nano-contact spin valve. Based on the 
derived analytical formula wc studied the power of the 
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FIG. 1. (a) A nano-contact spin valve is schematically shown, 
(b) A schematic illustration of the penetration of the mag- 
netic wall into the ferromagnetic electrode. The shaded circle 
represents the area of the nano-contact. 



current-induced oscillation and the thermal stability of 
the magnetic wall. 

The nano-contact spin valve we consider is schemati- 
cally shown in Fig. [T](a). To model the anti-parallel con- 
figuration of the spin valve, we applied fictitious magnetic 
fields of magnitude H to the top (bottom) electrode in 
the positive (negative) x-direction. The top and bottom 
electrodes were assumed to be square plates with width 
L, depth L and thickness d. The nano-contact was as- 
sumed to be a cylinder with radius a and height h. The 
electrodes and contact were assumed to be made of the 
same ferromagnetic material. The z-axis was taken to be 
parallel to the nano-contact and the electrodes parallel 
to the xy-plane. The origin of the coordinates was set to 
be the center of the nano-contact. 

The direction of the magnetization can be expressed 
by the polar angle 8 and the azimuthal angle (j>. The 
energy in the contact is given by 
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where A is the exchange stiffness constant. Because the 
exchange interaction energy is dominant in the nano- 
contactji 3 -, we neglected the demagnetization and Zeeman 
energies. 
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The energy of the top electrode is given by 
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where r = y/x 2 + y 2 , = <90/<9r, 9 = 89 /dr, M s is the 
saturation magnetization, and /io is the permeability of 
the vacuum. The first, second and third terms of Eq. ([2]) 
represent the exchange interaction energy, demagnetiza- 
tion energy, and Zeeman energy, respectively. The energy 
of the bottom electrode E-q is the same as Et- Assuming 
that the L is much larger than the penetration length of 
the magnetic wall, the system can be regarded to be ro- 
tationally symmetric about the z-axis; i. e., the angles 
and 9 depend only on the radius r. 

It is worth pointing out that, in the ground state, mag- 
netization in the top electrode m^r) and that in the 
bottom electrode m^{r) are related as (— mip, mij., m^) = 
(wj, ?7ig, trig) due to the symmetry of the system. In the 
following, we calculate the magnetic structure by taking 
this fact into account. 

Let us first consider the Bloch wall, where we assume 
that 9 = tt/2, d9/dz = 0, and 6 = 0. The energy of the 
contact is given by 
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To derive the last equation, we have assumed that 
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We hereafter call 



0o the boundary angle. We have also used the fact that 
is a linear function of z. Note that the corresponding 
Euler equation is d 2 <fi/dz 2 = 0. The energy of the top 
electrode is 
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The corresponding Euler equation was obtained as 



1 



2A 



0. 



(4) 



(5) 



The angle at r = a should be consistent with the 
boundary condition <t> z =h/2 = 0o for the contact and 
magnetic moments align in the direction of the magnetic 
field at r = oo. We therefore obtained the boundary con- 
ditions (j)\ r =a = 0o an d <p\ r =oo = 0. In the bottom elec- 
trode the boundary conditions should be <p\ r=a = ft — 0o 
and 4>(oo) = it. The energy of the electrode, which is a 
function of 0o , was obtained by substituting the solution 
of Eq. (O into Eq. Q and performing the integral. The 
boundary angle 4>q was obtained by minimizing the total 
energy E tota \ = E c + E T + E B . 

We introduce the characteristic length 2 ^ determined by 
the competition between the exchange interaction energy 
and the Zeeman energy as I = y / 2A/n Q M s H and use the 
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FIG. 2. (a) The energies of the contact E c , the top and 
bottom electrodes Et + Eb , and the total energy are plotted 
against the boundary angle (fro by the dotted, dashed, and 
solid lines, respectively, (b) The function G(a) of Eq. (|10[) 
is plotted as a function of 5. (c) The boundary angle 0o is 
plotted as a function of the thickness of the electrode d. (d) 
The analytical solution of the angle <j) i n the top electrode 
given by Eqs. ([8]) and is plotted by the solid line against 
the distance r from the center of the contact. The result 
obtained by the Monte Carlo simulation is plotted by the 
dotted line. The parameters are given in the main text. 



superscript to indicate the normalized values such 
as f = r/l, a = aft, d = d/l, and h = h/L Then Eq. © 
is expressed as 



1 • 



sin = 0, 
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where <j> = d 2 <fi/dr 2 and <p = d<p/df. Assuming that 
< 7r/4 in the top electrode, the last term of Eq. © can 
be approximated as sin ~ and we have 



^0-0 = 0. 

r 
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Equation ([7]) is the zeroth-ordcr modified Bessel equa- 
tion whose solutions are the zeroth-ordcr modified Bessel 
function of the first kind Io{r) and the second kind K (f). 
Since Io(r) diverges at the limit of f — > oo, the angle 0(r) 
should be expressed as 
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Substituting the approximation cos < 
Eq. © into Eq. Q we have 
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FIG. 3. (a) The twist angle of the Bloch wall is plotted as a 
function of the contact radius a and contact height h. The 
interval of the contours is 30° . (b) The same plot of the Neel 
wall. 



FIG. 4. (a) The difference between the twist angles 9* and <j>* 
is plotted against the contact radius a and the contact height 
h. (b) The same plot for the energy difference between the 
Neel and the Bloch walls. 



Here, K\{r) is the first modified Bessel function of the 
second kind. 

In Fig. [5] (a), wc plot the energies Eq (dotted-linc), 
Et + Eb (dashcd-line) and E to ta\ (solid-line) as func- 
tions of 4>o. Here and hereafter we use the following pa- 
rameters. The saturation magnetization is M s = 1.7 '4 x 
10 A/m, the exchange stiffness constant A = 20.7 pJ, 
the external field H = 7.96 x 10 4 A/m (1 kOe), the thick- 
ness of the electrodes d = 2 nm, the height of the contact 
h = 2 nm, and the radius of the contact a = 1 nm. The 
corresponding characteristic length is £ = 15.4 nm. Note 
that two parameters a and h will be varied in Figs. [2](c), 
l(a-b) andH(a-b). As shown in Fig. [1(a) E c {E T + E B ) 
is a parabolic function of the boundary angle 0o whose 
axis of symmetry is located at 0o = (tt/2), and -Etotai 
takes the minimum value at a certain value of 0o . Setting 
the first derivative of -Etotai to be zero, we found that the 
boundary angle of the ground state is 
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For the nano-contacts with a ~ h, one can easily see that 
the boundary angle 0o depends strongly on the value of 
d because the function G(a) ~ a as shown in Fig. [5] (b). 
In the limit of d — > oo the boundary angle 0o converges 
to zero, which means that the magnetic wall is perfectly 
confined in the contact as expected from the discussions 
in Ref. [H]. The value of 4>o increases with decreasing d 
and reaches ir/2 at the limit of d — > 0. In Fig. [2] (c) 
we plot the value of 0o as a function of the thickness of 
the electrode d. If we assume that the thickness of the 
electrode is the same as the height of the contact; i.e, 
d = 2 nm, the boundary angle is as large as 0o — 37.7°. 

We also performed Monte Carlo^ (MC) simulation to 
confirm the accuracy of the analytical formula we de- 
rived. In the MC simulation, the width and depth of the 
electrode L were set to be L = 100 nm, which is large 
enough to eliminate the size dependence. As shown in 
Fig. [5] (d) the angle 0(r) obtained by the MC simulation 
was well reproduced by our analytical formula of Eq. ([8]) 
with Eq. ([TTT) . 



Let us move on to the Neel wall where we assume that 
the azimuthal angle = 0, dtfi/dz = 0, and = 0. Wc 
also assume that — ir/2 < < tt/2. The energy of the 
contact is given by 
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6\ z =-h/2 = $0 — 7r/2. Note that Eq. (fl"2j) is equivalent to 
Eq. By rewriting in Eq. © with S6 = n/2-0 and 
taking up to the second order of 80, we obtain 
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where H c h = H + M$. The Euler equation for the Neel 
wall is given by 
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The boundary conditions are S0\ r=a = 0q and ^0| r _ ) . oo = 
in the top electrode, and S0\ r=a = n — and SOlr^oo = 
7T in the bottom electrode. Since Eq. (|14p takes the same 
form as Eq. ([7]) except that H is replaced with H c g, we 
can calculate SO in a similar way as we calculate the Bloch 
wall, by replacing the characteristic length £ with £ c ff = 
^2A/ LiQAI s H c ff . Assuming the same parameters as in 
Fig. [2] (a) £ s is estimated to be 3.1 nm. The energy 
of the top electrode Et and the boundary angle 6$ of 
the ground state of the Neel wall are calculated from 
Eqs. ((9j) and (|11[) . respectively, by changing the length 
for the normalization from £ to £ g. 

In order to discuss the difference of the magnetic struc- 
tures between the Bloch wall and the Neel wall we intro- 
duce the twist angles 4>* and 0*, which are defined as 
<t>* = \4>\z=h/2 - 0U=-va| = 7T - 20o for the Bloch wall 



and 0* ee 



\z=h/2 " 



\z=-h/2 



= n - 20 for the Neel wall. 



As shown in Figs. [3] (a) and (b) 0* is always larger than 
0* because the penetration of the Neel wall costs more 
demagnetization energy in the electrodes compared with 
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that of the Bloch wall. Since the resistance of the mag- 
netic wall is proportional to the twist angled, the con- 
ditions for maximizing the power of the current-induced 
oscillation between the Bloch wall and the Neel wall are 
the same as those for the difference between 0* and 4>*. 
In Fig. |4] (a) we plot the difference 8* — (j>* as a function 
of a and h. One can see that the optimum conditions 
for maximizing the power of the STO based on the nano- 
contact spin valve are approximately given by a ~ h. 

Figure [4] (b) shows the difference between the en- 
ergies of the Neel wall and the Bloch wall defined as 
AE = -EnccI — -E-Bioch as a function of a and h. The en- 
ergy difference AE, and therefore the thermal stability of 
the Bloch wall, increases with increasing a or decreasing 
h. In the plotted region (0 < a, h < 5nm) the Bloch wall 
is the ground state and the Neel wall the excited state, 
contrary to the results of the perfectly confined magnetic 
wall shown in Ref. [T^. The energy difference AE for the 
perfectly confined magnetic wall originates from the de- 
magnetization energy in the contact that we neglected. 
Following Ref. [TH we can estimate AE ~0.16 cV for a 
perfectly confined magnetic wall with a =1 nm and h =2 
nm. On the other hand, for the the nano-contact spin 
valve with thin ferromagnetic electrodes, the main con- 
tribution to AE is from the demagnetization energy in 



the electrodes. The estimated value of AE is about 3.0 
eV, which is about 19 times larger than that of the per- 
fectly confined magnetic wall. 

In summary, we theoretically studied the effects on a 
magnetic structure of a nano-contact spin valve with thin 
ferromagnetic electrodes. We derived analytical formulas 
for the magnetic configurations of the Bloch wall and the 
Neel wall, and the corresponding energies. We showed 
that, compared with the Bloch wall, the penetration of 
the Neel wall into the electrodes is suppressed by the 
increase of the demagnetization energy. We found the 
optimum conditions for maximizing the power of the STO 
based on the nano-contact spin valve are given by a ~ h. 
We also found that the thermal stability of the magnetic 
wall increases as the nano-contact's radius increases or 
height decreases. 
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